Pseudogap phase of high-Tc compounds described within the LDA+DMFT+2 approach 
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LDA+DMFT+Sk approach was applied to describe pseudogap phase of several prototype high-Tc compounds e.g. hole doped 
Bi2Sr2CaCu208-^ (Bi2212) and La2-xSr;,Cu04 (LSCO) systems and electron doped Nd2-xCe;cCu04 (NCCO) and Pr2-xCe;cCu04 
(PCCO), demonstrating qualitative difference of the Fermi surfaces (FS) for these systems. Namely for Bi2212 and LSCO the so 
called "hot-spots" (intersection of a bare FS and AFM Brillouin zone (BZ) boundary), where scattering on pseudogap fluctuations is 
most intensive were not observed. Instead here we have Fermi arcs with smeared FS close to the BZ boundary. However for NCCO 
and PCCO "hot-spots" are clearly visible. This qualitative difl'erence is shown to have material specific origin. Good agreement 
with known ARPES data was demonstrated not only for FS maps but also for spectral function maps (quasiparticle bands including 
lifetime and interaction broadening). 
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1. Introduction 

One of the most prominent phenomena in high-Tc cuprates 
physics is the so called pseudogap |1]. Here we present 
an overview of our recent works Refs. 10, S B S] on 
LDA-hDMFT-hEk computational scheme applcations. This 
I scheme is generalization of dynamical mean-field theory 
DMFT [6] and LDA-hDMFT |7] (LDA - local density approx- 
imation) approach allowing to include non-local scale depen- 
dent eff'ects. To include pseudogap fluctuations efl'ects 
important for cuprate physics we supplied (in additive man- 
ner) conventional DMFT with an external k-dependent self- 
energy Ek. For the pseudogap state describes the interac- 
tion of correlated electrons with non-local (quasi) static short- 
ranged collective Heisenberg-like AFM or SDW-like spin fluc- 
tuations (id,^. 

Within LDA-hDMFT-hSk approach several high-T^ prototype 
jcompounds e.g. hole doped Bi2Sr2CaCu208-^ (Bi2212) |2] 
and La2-xSr;cCu04 (LSCO) |3] as well as electron doped 
Nd2-xCe;,Cu04 (NCCO) il and Pr2-;cCe;,Cu04 (PCCO) |5] 
were studied. Since most powerful experimental tool to ac- 
cess electronic properties of the pseudogap state is angular re- 
solved photoemission spectroscopy (ARPES) |12, 13, 14] we 
performed comparison of LDA-hDMFT-hEk calculated spectral 
functions and Fermi surfaces with available ARPES quasiparti- 
cle bands and Fermi surface maps. Two-particle properties can 
also be described by this approach ifisl] . e.g. calculated optical 
spectra in the pseudogap state compare well with experimental 
data for Bi2212 L2J and NCCO 14J. 



2. LDA+DMFT+i: computational details 

Crystal structure of Bi2212 ||2|], NCCO H and PCCO H 
has tetragonal symmetry with the space group I4/mmm, while 
LSCO has ortorhombically distorted structure Bmab |3]. For 
further crystallographic data used within our LDA-hDMFT-hEk 
approach see Refs. 10, S H 0]- Well known quasi two- 
dimensional nature of these compounds determines its physi- 
cal properties. Physically most interesting are the Cu02 layers. 
Those layers provide antibonding Cu-3d(x^-y^) partially filled 
orbital, whose dispersion crosses the Fermi level. Thus we are 
using this eff'ective LDA calculated Cu-3d(x^-y^) antibond- 
ing band as a "bare" band in LDA-^DMET-hSk computations. 
Corresponding hopping integral values obtained within the lin- 
earized muflftn-tin orbitals (LMTO) method 1 16] and further ap- 
plication of the A^-th order LMTO (NMTO) approach tlTJ are 
listed in Table 1 . 

Next to perform DMFT calculations one should set up on-site 
Coulomb interaction values. The values of Coulomb interaction 
on eff'ective Cu-3d(x^ -y^) orbital U obtained via constrained 
LDA computations! 18] are also presented in the Table 1. 

To account for the AFM spin fluctuations, a two-dimensional 
model of the pseudogap state is applied. [HI] Correspond- 
ing k-dependent self-energy Sj^ ill [13, [HI] describes nonlo- 
cal correlations induced by (quasi) static short-range collective 



Table 1 : Calculated energetic model parameters (eV) and experimental correla- 
tion length ^. First four Cu-Cu in plain hopping integrals t, f , f , t'" , interplain 
hopping value t^, local Coulomb interaction U and pseudogap potential A. 
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Figure 1: LDA+DMFT+Ek spectral functions for Bi2212 (upper panel) and 
NCCO (lower panel) along of noninteracting FS in 1/8 of BZ. Dashed-black 
line corresponds to "hot-spot" (Ref. | 4]) 



Heisenberg-like AFM spin fluctuations. 11 19l] 

The Sj^ definition contains two important parameters: the 
pseudogap energy scale (ampHtude) A, representing the energy 
scale of fluctuating SDW, and the spatial correlation length ^. 
The latter is usually determined from experiment. The A value 
was calculated as described in Refs. The value of cor- 

relation length was taken in accordance with the typical value 
obtained in neutron scattering experiments on NCCO | 21] and 
LSCO 1 22]. Employed values of A and ^ for all considered 
systems are shown in Table 1 . To solve DMFT equations nu- 
merical renormalization group (NRG, Refs. EQ) was em- 
ployed as an "impurity solver". Corresponding temperature of 
DMFT(NRG) computations was 0.01 1 eV and hole or electron 
concentrations were 15%. 



3. Results and discussion 

Based on extended analysis of LDA-hDMFT-hSk results and 
experimental ARPES data the origin of pronounced "hot-spots" 
(cross-point of the Fermi surface and umklapp surface) for elec- 
tron doped systems was established. Also it was shown 
that hole doped systems have only Fermi arcs Ei. Fig.m 
displays LDA-FDMFT-hSk spectral functions along 1/8 of non- 
interacting FS from the nodal point (top curve) to the antin- 
odal one (bottom curve). Data for Bi2212 is given in left panel, 
NCCO — right panel of Fig. [T] For both compounds antin- 
odal quasiparticles are well-defined — sharp peak close to the 
Fermi level. Going to the nodal point quasiparticle damping 
grows and peak shifts to higher binding energies. This behav- 
ior is confirmed by experiments Refs. |25, 26] (for comparison 
with experiment see Ref. 0]). Let us interpret the spectral func- 
tion peaks based on the LDA-FDMFT-hSk results. Namely, for 
Bi2212 nodal quasiparticles are formed by low energy edge of 
pseudogap, while for NCCO they are formed by higher energy 
pseudogap edge. Also in NCCO there is obviously no bilayer 
splitting eff'ects seen for Bi2212 (left panel of Fig. [T]). 

"Hot-spots" for NCCO are closer to the BZ center ||4|]. In 
Fig. [T] one can see it from the position of the dashed-black line 



which corresponds to the "hot-spot" k-point. For Bi2212 scat- 
tering from neighboring BZ amplify each other and instead of 
just"hot-spot" we see rather extended "destructed" Fermi sur- 
face towards the BZ boundaries. Such strong scattering comes 
from scattering processes with momentum transfer of the order 
ofQ=(7r,7r)UriB lin . correspondmg to AFM pseudogap fluc- 
tuations. Qualitatively the same picture is found also in LSCO 
(see Fig. [21). 

Recent experimental and theoretical LDA-hDMFT-hSk Fermi 
surface maps | 3] are shown in Fig. Oat panels (a) and (b) corre- 
spondingly. Both pictures reveal strong scattering around (7r,0)- 
point which we associate with scattering in the vicinity of the 
so-called "hot-spots" which are close to the (7r,0) 10,13]. Along 
nodal directions we observe typical Fermi arcs. They are pretty 
well seen in the theoretical data while in experiment we ob- 
sereve just narrow traces of them (Bi2212 Fermi surface is com- 
pared with experiment in Ref. |2]). 

Another possibility to compare LDA-l-DMFT-l-Sk results 
with ARPES data is spectral function colour maps plotted 
along symmetry lines. In Fig. [3] we present LDA-hDMFT-hEk 
intensity plots along the high symmetry lines for NCCO 
(upper panel) in comparison with high-energy bulk sensi- 
tive angle-resolved photoemission data of Ndi.85Ceo.i5Cu04 
(lower panel). 1 4] Indeed we see quite a good agreement of 
LDA-hDMFT-hSk and experimental data. For the M - F direc- 
tion there is not very much going on. Basically we see both in 
theory and experiment very intensive quasiparticle band. For 
the M-T direction less intensive shadow band is not resolved 
in the experiment. 

More interesting situation is observed for T-X-M directions. 
At F-point there is band in the experiment starting at about -1.2 
eV. It is rather intensive and goes up in energy. Suddenly there 
is almost zero intensity at about -0.3 eV. Then in the vicinity 
of the X-point intensity rises up again. In the X- M direc- 
tion around -0.3 eV on the right side of X-point there is also 
quite intensive region. At a first glance one can think that it is 
the same band with matrix element eff'ects governing intensity. 
However based on analysis of Ref. |4] one can conclude that 
this low intensity region is the forbidden gap between shadow 
and quasiparticle bands. The "horseshoe" around X-point is 
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Figure 2: Fermi surfaces of LSCO at x=0.14 from experiment (left panel) 
and LDA-i- DMFT+Ek computations (right panel) Red crosses on the left panel 
correspond to experimental kf values. (Ref. [SJ) 
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Figure 3: Comparison of LDA+DMFT+Ek spectral functions (upper panel) 
for NCCO along BZ high- symmetry directions with experimental ARPES 
(Ref. U) (lower panel). 



formed by the shadow band on the left and the quasiparticle 
band on the right for upper branch and other way round for 
the lower branch. As a consequence of that there is also inten- 
sive shadow FS sheets around (n/afl) point. Rather intensive 
nondispersive states at about -1.0 eV within experimental data 
can be presumably associated with the lower Hubbard band and 
possible admixture of some oxygen states. Let us also suppose 
that high intensity at -0.3 eV for X point may be interpreted 
not as a van-Hove singularity of bare dispersion but rather of 
high-energy pseudogap branch |4]. 

One more fascinating comparison for LDA-hDMFT-i-Ej^ re- 
sults with experimental ARPES data is recently reported by 
us in the Ref. ||5]. In Fig. |4] an extended picture of PCCO 
Fermi surfaces is presented (panel (a) — LDA-hDMFT-hSj^ re- 
sults, panel (b) — experimental ARPES data). Strictly speak- 
ing Fig. |4] is a color map in reciprocal space of the corre- 
sponding spectral function plotted at the Fermi level. FS is 
clearly visible as reminiscence of non-interacting band close 
to the first Brillouin zone border and around (7r/2,7r/2) point 
(so called Fermi arc), where the quasiparticle band crosses the 
Fermi level. There is an interesting physical eff'ect of partial 
"destruction" of the FS observed in the "hot-spots", points that 
are located at the intersection of the FS and its AFM umklapp 
replica. This FS "destruction" occurs because of the strong 
electron scattering on the antiferromagnetic (AFM) spin (pseu- 
dogap) fluctuations on the copper atoms. Also the "shadow" 
FS is visible as it should be for AFM folding. As no long-range 
order is present in the underdoped phase the "shadow" FS has 
weaker intensity with respect to FS. The PCCO FS is very sim- 
ilar to that of Nd2-xCe;cCu04 (NCCO), which belongs to the 
same family of supperconductors ||4l |25|] . 

Let us compare theoretical (upper panels) and experimental 
(lower panels) energy quasiparticle dispersion for most charac- 
teristic cuts introduced in Fig. [4] (see Fig. O. Theoretical data 
were multiplied by the Fermi function at a temperature of 3 OK 




Figure 4: (a) Extended Fermi surfaces for PCCO — LDA+DMFT+Ej^ 
data. White rectangle on panel (a) schematically shows the part of recipro- 
cal space measured experimentally (panel b). Lower left corner is X-point 
(;r,0).(Ref. |5J) 

and convoluted with a Gaussian to simulate the efl'ects of exper- 
imental resolution, with further artificial noise added. 
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Figure 5: Energy-momentum intensity distributions for the specific cuts drawn 
in Fig. |4] (upper panels — theoretical data, lower panels — experimental pho- 
toemission intensity). To judge about the absolute intensities of the "shadow" 
(1) and main band (2) cut 1 contains an MDC curve integrated in an energy 
window 60 meV centered at the Fermi level (FL). Similarly integral EDC for 
cut 2 ("hot-spot") shows suppression of the intensity at the FL as compared to 
cut 3, which is located further away from the "hot-spot". (Ref. |4]) The FL is 
zero. 

The Cut 1 intersects quasiparticle and "shadow" Fermi sur- 
faces close to the Brillouin zone border. One can find here 
a "fork"-like structure formed by the damped "shadow" band 
(-0.5-0 arb.u.) and better defined quasiparticle band (0.5-1 
arb.u.). This structure corresponds to preformation of FS cylin- 
der around (7r,0) point. The Cut 2 goes exactly through the 
"hot-spot". Here we see a strong suppression of the quasipar- 
ticle band around the Fermi level similar to NCCO as shown 
in Fig. O The Cut 3 crosses the Fermi arc, where we can see 
a very well defined quasiparticle band. However weak inten- 
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sity "shadow" band is also present. For the case of long range 
AFM order and complete folding of electronic structure, FS and 
its "shadow" should form a closed FS sheet around (7r/2, 7r/2) 
point, while in the current case the part of the pocket formed 
by the "shadow" band is not so well defined in momentum 
space. As can be seen there is a good correspondence between 
the calculated and experimental data in terms of the above de- 
scribed behavior, which is also similar to the results reported 
for Nd2-xCe;,Cu04 (NCCO) in our earUer work.fl 

4. Conclusion 

Here we summarize our recent results on LDA+DMFT+Sj^ 
investigations of pseudogap state for a number of copper high- 
Tc compounds. We considered for the main prototype sys- 
tems: hole doped - Bi2212 [M]^ and LSCO Mi; electron 
doped - PCCO |5] and NCCO [4]. For all compounds the 
LDA-hDMFT-i-S]^ calculations show that Fermi-liquid behav- 
ior is still conserved far away from the "hot-spots" (antinodal 
direction), while the destruction of the Fermi surface observed 
in the vicinity of "hot-spots" (close to nodal direction). This 
destruction is due to strong scattering of correlated electrons on 
short-range antiferromagnetic (pseudogap) fluctuations. More- 
over the origin of clearly observed "hot-spots" for electron 
doped systems (in contrast to hole doped ones with a Fermi 
arcs only) is established. Comparison between experimental 
ARPES and LDA-hDMFT-hSj^ data reveals a good semiquanti- 
tative agreement. The experimental and theoretical results ob- 
tained once again support the AFM scenario of pseudogap for- 
mation not only in hole doped HTSC systems |,2fci3J but also in 
electron doped ones . 
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